Topological Fermi arcs in superfluid 3 He 
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We consider fermionic states bound on domain walls in a Weyl superfluid 3 He- A and on interfaces 
between 3 He-A and a fully gapped topological superfluid 3 He-B. We demonstrate that in both 
cases fermionic spectrum contains Fermi arcs which are continuous nodal lines of energy spectrum 
terminating at the projections of two Weyl points to the plane of surface states in momentum space. 
The number of Fermi arcs is determined by the index theorem which relates bulk values of topological 
invariant to the number of zero energy surface states. The index theorem is consistent with an exact 
spectrum of Bogolubov- de Gennes equation obtained numerically meanwhile the quasiclassical 
approximation fails to reproduce the correct number of zero modes. Thus we demonstrate that 
topology describes the properties of exact spectrum beyond quasiclassical approximation. 

PACS numbers: 



I. INTRODUCTION 

Chiral Weyl fermions represent the fermionic sector 
in the Standard Model of particle physics. The point 
nodes in the spectrum of chiral quarks and leptons are 
topologically protected, as well as their condensed mat- 
ter counterparts, which are called Dirac or Weyl points 
(on topology of Weyl points in particle physics and con- 
densed matter see e.g. Refs. Q]-|5|). 

The nodal topological systems with Weyl fermions 
demonstrate different types of the bulk-surface and the 
bulk-vortex correspondence. Due to the bulk-vortex cor- 
respondence, the cores of some vortices in the Weyl su- 
perfluid 3 He-A contain a dispersionless branch of bound 
states with zero energy which is one-dimensional flat 
band, discussed first by Kopnin and Salomaa in 1991 6 . 
The end points of this flat band arc determined by the 
projections of Weyl points to the direction of vortex axis 
(see Refs. jjP for the topological origin of this fiat band, 
and Refs. l9Ml2lfor discussion of the topological fiat bands 
in general). Due to the bulk-surface correspondence, the 
surface of a system with Weyl points contains another ex- 
otic object - the Fermi arc - the ID Fermi line in the 2D 
momentum space, which terminates on the projections of 
two Weyl points to the plane of the surface. The Fermi 
arc on the surface of 3 He- A shown in FigJ5^ has been 
considered in Ref. and that on the surface of topolog- 
ical semi-metals with Weyl points - in Refs. Il4rjl6l The 
flat band in the vortex core and the Fermi arc on the sur- 
face are momentum-space analogs of the real-space Dirac 
string terminating on two magnetic monopoles 9 . 

Fermi arcs appear also at the interface, which separates 
two degenerate states of 3 He-A with the opposite direc- 
tions of the orbital anisotropy vector 1 (in chiral super- 
fluid 3 He-A vector 1 determines the orientation of a spon- 
taneous orbital angular momentum of this anisotropic 
liquid and also determines the positions p = ±ppl of 
two Weyl points on the Fermi surface with topological 
charges N3 = ±2, if spin degeneracy is taken into ac- 



count). Bound states emerging at one of the representa- 
tives of such interface - the continuous 1-soliton - have 
been calculated in Ref. [l?]- Here we consider bound 
states emerging on the singular domain wall in 3 He-A 
discussed in Refs. GHH Zero energy edge states on such 
domain wall in a thin film of 3 He- A and the topological 
bulk-edge correspondence for this 2+1 system were con- 
sidered in Ref. [2CJ. The Fermi arc emerging in the 3+1 
system forms a special configuration in momentum space, 
see Fig. [2p. We also consider the Fermi arc emerging at 
the interface between two topologically different quantum 
vacua: the 3 He-A with Weyl points and the fully gapped 
3 Hc-B, which also has nontrivial topology in momentum 
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(AB interface). 



space 

In the present paper we consider the Fermi arcs in su- 
perfluid 3 He solitons which differ qualitatively from that 
existing in other systems such as the surface of 3 He-A and 
Weyl semimetals. The reason is that they can be ob- 
tained ultimately beyond the quasiclassical approxima- 
tion. It was found that the quasiclassical approximation 
yields a large number of spectrum branches which inter- 
sect the Fermi leveh 17 ' 25 The number of such branches 
depends on the parameters of the solition in contradic- 
tion with the topological index theorem. We show that 
the reason for this contradiction is that the quasiclassical 
Bogoluibov - de Gennes (BdG) equations inherently miss 
the normal reflection of quasiparticles from the spatially 
inhomogeneous superfluid order parameter and take into 
account only Andreev reflection. Thus the system of qua- 
siclassical BdG equations is often called the system of 
Andreev equations. As we will discuss below the nor- 
mal reflection is crucial to describe the Fermi arcs in 
superfluid 3 He solitons. In this case fermionic zero en- 
ergy states are realized on the quasiparticle trajectories 
passing at the small sliding angle to the soliton plane. 
For such trajectories the normal reflection is strongly en- 
hanced and change qualitatively the behavior of spectral 
branches near the zero energy points on Fermi arcs. Here 
we find numerically exact solutions of the spectral prob- 
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lems for BdG equations which confirms the topological 
index theorem predictions. Thus our results demonstrate 
that topological bulk-edge correspondence describe the 
properties of the exact spectrum of localized boundary 
states which in some cases are missed in the quasiclassi- 
cal approximation. 

The structure of the paper is as follows. In section 
ITT1 we introduce the model which is the BdG equations 
and the order parameter structures of domain wall in A 
phase and AB interface. The results are presented in sec- 
tion Hn] where we consider at first topological nature of 
Fermi arcs and then the calculations of fcrmionic spec- 
trum. We discuss the failure of quasiclassical approxima- 
tion to describe Fermi arcs in the fermionic spectrum of 
domain wall in A phase and AB interface in section HVl 
We summarize the results in section IVl 



II. THE MODEL 

We calculate the spectrum of eigenstates of Bogolubov- 
de Gennes (BdG) Hamiltonian describing spin-triplet p- 
wave superfluids/superconductors 

H = [e(p) - /i]f 3 + fiReA - f 2 ImA, (1) 

where e(p) = (p x + p y +pl)/2m, /i is chemical poten- 
tial, pi = —iVi and n are Pauli matrices of Bogolubov- 
Nambu spin acting on the wave function ip = (u, v) T with 
particle and hole components u and v. The gap operator 
is A = A a iU a pi/pF, where a a are Pauli matrices of spin 
(in 3 He it is nuclear spin), pp — y/2mjl is the Fermi mo- 
mentum. We will consider the spectral problem of BdG 
Hamiltonian (JTJ) for the fermionic states bound on the ID 
solitons in superfluid 3 He namely the domain wall of the 
A phase and the AB interface. 

The order parameter in superfluid 3 He is 3 x 3 ma- 
trix A ai where the Greek and Latin indices correspond 
to the spin and orbital variables. We will consider the 
fermionic modes localized on different ID order parame- 
ter solitons, namely the domain wall of a 3 He-A and the 
AB interface with order parameter inhomogeneity along 
the vector n = x so that the momentum projections p y 
and p z are conserved. At first we introduce the order 
parameter structure corresponding to the discussed soli- 
tons. 

Domain wall of A phase. We consider the configu- 
ration of one of the domain walls possible in 3 He-A ac- 
cording to the symmetry classification^- the one which 
is shown in Fig. QJi. The orbital vector 1 of the order 
parameter points down (1 || — z) as x —} — oo and points 
up (1 || z) as x — > +oo. This domain wall configuration 
can be approximated by the following Ansatz: 

A ai (x) = A d a [x l + if(x)yi] , (2) 

where the unit vector d represents the spin part of the or- 
der parameter, which is fixed if the spin-orbit interaction 
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FIG. 1: (Color online) Model spatial dependencies of the or- 
der parameter components corresponding to the most sym- 
metric and energetically preferable ID t opolog ical defects in 
superfluids 3 He according to Refs.(18 19.31 32): (a) domain 
wall of the A phase and (b) AB interface. The order param- 
eter is normalized to the bulk value. 

is neglected, and we put d = x. Then the gap function 
is given by 

A A a = c a x [p x + if{x)p y ] , (3) 

where cq — Ag/pp is the parameter which plays the 
role of the longitudinal speed of light in bulk 3 He- A 2 , 
f(x) is arbitrary monotonic function with asymptotic 
/(±oo) = ±1. The topology of bound states does not de- 
pend on the details of the function f(x) in Eq.([3]) if it is 
monotonic. Note that the same model (modulo unimpor- 
tant d vector orientation) describes domain wall between 
chiral domains in p + ip superconductors like Sr2RuOi 2i 
thus the fermionic spectra are identical in these cases. 

AB interface. Let us fix the order parameter on the 
B phase side of the interface, i.e. at x = +oo, in the 
following form 

A ai (x = +oo) = A B (x a Xi + y a yi + z a Zi). (4) 

The configuration of the A phase at x — > — oo can be 
different for different realizations of the domain wall (see 
e.g. Ref. 31). They are described by the relative ori- 
entations of the vectors 1, d and n (the normal to the 
AB interface which in our case is n = x). First we shall 
consider the case d = x, 1 = z which has the lowest 
energy 3 ^ 

A aj (x = -oo) = A A x a {xj ~ iyj). (5) 

We model the domain wall at the AB interface with A^ = 
Ab = A and the switching between bulk phases (|4I5[) 
as follows 

A xx = A Q = const (6) 
A yy = A zz = Ao/i(x) (7) 
A xy = -iA f 2 (x) (8) 

where fi^fa) are arbitrary monotonic functions with 
asymptotics f\{— oo) = 0, /i(+oo) = 1 and f%(— oo) = 1, 



3 



/2(+oo) = 0. The order parameter components ([6]) are 
shown in Fig.(QJ>). Then the gap operator is given by 



A 



AB 



cqCTx [Px - ifi{x)Py] + c fi{x) [ayPy + a z p 



III. RESULTS 
A. Topology of Fermi arcs in superfiuid 3 He 
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FIG. 2: (Color online) The manifold of zero energy states in 
Py,Pz plane in the spectrum of bound fermions e(p y ,p z ) = 
forms (a) the Fermi arc (solid blue line) on the boundary 
of 3 He-A according to Ref. [l3|; (b) two Fermi arcs on the 
domain wall in Eq.Q. Only single spin projection is con- 
sidered. Thick arrows show directions of Fermi velocity at 
the Fermi arcs. In (a) the Fermi arc has topological charge 
N = +1, which satisfies the index theorem following from the 
bulk-surface correspondence and the momentum space topol- 
ogy of Weyl points in bulk 3 He-A. Fermi arc terminates on 
the projections of Weyl points to the surface of 3 He-A. The 
spectrum of bound states terminates at the dashed line where 
the spectrum merges with the bulk spectrum. The region of 
continuous spectrum is shown by yellow shading. In (b) at 
Pz = the Fermi velocity at the Fermi arcs is in the same 
direction, vf = —VfY which demonstrates that both Fermi 
arcs have the same topological charge N = +1, which to- 
gether satisfy the index theorem N$(left) — Nz{right) — +2, 
in agreement with momentum space topology of Weyl points 
in bulk 3 He-A on two sides of the domain wall in Fig. [3] 
Fermi arcs terminate on the projections of Weyl points to the 
interface. The spectrum of bound states has discontinuity at 
p y = where the spectrum merges with the bulk spectrum 
(the edge of continuum is shown by dashed line). 



FIG. 3: Topology of the fermionic bound states on the domain 
wall in 3 He-A. Momentum space topology of Weyl points in 
bulk 3 He-A on two sides of the wall prescribes existence of 
the Fermi arcs in the spectrum of the fermionic states in the 
soliton or at the interface between the bulk states with dif- 
ferent positions of Weyl points. In the considered case of the 
domain wall the Weyl points on two sides of the interface 
have the same positions in momentum space, but opposite 
topological invariants N3. This leads to two Fermi arcs ter- 
minating on the projections of the Weyl points to the plane 
of domain wall (Fig. [2b) according to the index theorem 
-/V3 (left) — N3 (right) = 2. This is distinct from the single 
Fermi arc on the surface of the 3 He-A in Fig. [2Jl. 



[2^,. Topological origin of two Fermi arcs on the domain 
wall is demonstrated in Fig. |3] The Weyl points on two 
sides of the interface have the same positions in momen- 
tum space, but opposite values of topological invariant 
N 3 . 

The bulk spectrum in the plane with fixed p z in mo- 
mentum space has no nodes if \p z | + 1 pf and thus corre- 
sponds to the spectrum of 2D insulator. This insulator is 
topological for \p z \ < pp, since it is described by nonzero 
topological invariant introduced in Refs. [27H3G1 for 2+1 
systems: 



N 3 (p z ) 
1 

~ 4^2 



tr 



dp x dp y dw Gd P:r G~ 1 Gdp v G~ 1 Gd UJ G~ 1 



(10) 



As we have mentioned in the Introduction the spec- 
trum of fermions bound to the surface of superfiuid 3 He 
A is known to contain the Fermi aroi 3 -. Let us now com- 
pare the Fermi arcs on a domain wall in 3 He-A with the 
Fermi arc on a surface of 3 He- A shown in Figs. ((2^,) and 
(J2b) correspondingly. In both cases Fermi arcs originate 
from the Weyl points in the bulk spectrum. But in case 
of the domain wall the Weyl points exist in the bulk liq- 
uids on both sides of the wall, as a result there are two 
Fermi arcs in Fig. (2b instead of a single Fermi arc in Fig. 



Here G is the Green's function matrix, which in our 

non-interacting models is G = (iui — H^j . One has 

-^3(1^2 1 < pf) = +1 on one side of the wall and 
•^3(1^1 < Pf) = —1 on the other side. According to 
the index theorem^^, the difference between these two 
values determines the number of the zero modes at the 
interface between the 2+1 topological insulators for each 
\Pz\ < Pf- As a result one has two Fermi arcs at the 
soliton wall. 
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B. Calculation of the spectrum 

At first we note that the spectral problem for the do- 
main wall can be significantly simplified since the BdG 
Hamiltonian (fTJ) with the gap operator given by Eq.Q 
commutes with fermionic spin a x . Thus we transform 
the quasiparticle wave function tp = (u, v) T to remove 
the spin dependence of the order parameter u = a x u, 
v = v which yields the gap operator in the form 

A = Co \p x + if(x)p y ] . (11) 

Furthermore the spectral problem for the AB interface 
with the gap operator given by Eq.® can be mapped 
on the case of domain wall of the A phase. Let us 
transform the quasiparticle wave function components 
as follows u — a x u, v — v. We choose the spin basis 
Xa = (iPz,c±P± — Py) of the eigenstates of the opera- 
tor a± = (T z p y — o y p z corresponding to the eigenvalues 

<7j_ = ±1 where p± — \jp y + p\ ■ Then the order param- 
eter is diagonal in spin space and has the form of Eq.(fTT|) 
where 

f(x) = <rxpj_fi(x) -p y f 2 {x). (12) 

In order to obtain the bound states the function (fl2|) 
has to satisfy the condition /(+oo)/(— oo) < which 
yields cr±p y > 0. Thus in contrast to the case of the A 
phase domain wall considered above the spin degeneracy 
is removed since the proper spin state is determined by 
the condition of bound state existence. 

The continuous part of the spectrum e c (p) of fermionic 
excitations is determined by the eigen states of Hamilto- 
nian ([1]) at the bulk regions \x\ 3> r 

£c(p) = ±^[e( P )-^+cl(pl+pl). (13) 

The states localized at the domain wall are characterized 
by a single discrete quantum number which enumerates 
the energy branches and by the two continuous quan- 
tum numbers which are the projections of the quasipar- 
ticle momentum p y ^ z onto the domain wall plane. The 
energy of localized states is confined within the region 
|e| < min Px e c (p). 

In general the eigenvalue problem of BdG Hamilto- 
nian (fTJ) yields a system of differential equations which 
can not be solved analytically. However there are ap- 
proximate methods which can help to study qualitative 
features of the spectrum. First of all we will employ the 
semiclassical approach 2 ^ when the momentum operator 
is approximated by a number p x = p x {x) thus turning 
the differential BdG equations into the algebraic ones. In 
this case the spectrum of bound states is determined by 
Bohr-Sommerfeld quantization of classical periodic mo- 
tion between reflection points discussed in the Section 
IHIB II The semiclassical approximation fails when the 
distance between reflection points is too small. Particu- 
larly important is the case when two Andreev reflection 



points come close together forming the bound state which 
correspond to the so called 'zeroth branch' of the energy 
spectru m 17 i 25 . This case can be treated with the help 
of another approach - the quasiclassical approximation 
when one can use reduced order Andreev equations for 
the envelope wave functions discussed in Section MI B 21 
The quasiclassical approximation in BdG equations is ap- 
plicable as long as the normal reflection can be neglected. 
For the structures considered in the present paper the 
normal reflection is crucial to describe properly the zero 
modes of the fermionic spectrum. This can be done only 
with the help of numerical solution of the exact BdG 
system of equations discussed in Sectio nTlI B 31 

Hereafter we will choose the model form of Eq.(fTTj) 
with f(x) = tanh(a;/r) where r ~ £ is width of the 
domain wall and £ = hvp/^-o is the coherence length 
where vp — Pf/ 1 ^ 1 is the Fermi velocity. The fermionic 
spectrum in this configuration was considered in Ref. [25| 
in quasiclassical approximation. Here we will implement 
an analysis of the fermionic spectrum beyond the quasi- 
classics in order to study the the Fermi arcs supporting 
Majorana states. 

1. Semiclassical spectrum. 

We begin the analysis of the spectrum of localized 
states from the semiclassical approximation. Provided 
the condition \p x \^ 3> 1 is valid we use the expression for 
the energy (IT31 substituting the local value of the order 
parameter. In this case we obtain 

e = ±^J[pl/2m - ^ + c^pl + p(x)pl}. (14) 

where p, x = p,— (py+pl)/2m. From Eq. (fT4|) the function 
Px = Px{x) can be found which describes semiclassical 
orbits in (p x ,x) phase space. In general the orbits can 
have two types of stationary points determined by the 
nature of quasiparticle reflection. 

The normal reflection occurs at the points x = x™ 2 
defined by p x (xi 2 ) — or equivalently 

eoft,/« 2 ) = ±yje 2 -e 2 n 2 (15) 

and exists at the energy interval e n \ > \e\ > s n 2 where 

e„i = V ( c oPy) 2 + an d £ n2 = \px\- The Andreev re- 
flection occurs at x — x I 2 where Im(p x ) becomes non- 
zero 

C0Py/(s? j2 ) = ±y/e 2 -el 2 (16) 

and exists if e a i > |e| > £ a 2 where e a i = 
°°\JPy + 2m A t a; - (mco) 2 and e a2 = c oy /2mp, x - (mc ) 2 . 

One can see that e n \ > e a \ and e n2 > £ a 2 therefore 
Eq. (ll4[) determines two qualitatively different types of 
the enclosed quasiparticle orbits in (p x , x) space. That 
is for the energies (i) e n \ > \e\ > e a i the orbits have 



5 



only normal reflection points and for (ii) e n i > \e\ > £ a 2 
only Andreev reflection points. The orbits of type (i) and 
(ii) are shown in the Fig.(j4ji) by red dashed and green 
solid lines correspondingly. Provided (iii) e a \ > e n 2 there 
is also the third regime when the orbit has both normal 
and Andreev reflection points which is shown by the blue 
dash-dotted line in the Fig. (01. The energy spectrum is 
determined by the Bohr-Sommerfeld quantization 

j> p x dx = 2?r(n s + 7) (17) 

where n s is integer, 7=1/2 when the number of turning 
points is ever* 2 ^. Semiclassical spectral branches for n s — 
1,2 are shown in Fig.fQJ}) by solid lines. 

From the above semiclassical consideration we obtain 
that the energy of localized states is not bounded as func- 
tion of p y and p z . Indeed let us consider the case (i) and 
assume y 'e 1 — fi x -C co[p v | -C |e|. In this case we can 
put f(x) = x/r in the Eq. (fT4|) which yields then the 
harmonic oscillator spectrum 

E ns (p y ,p z ) = \fi x \ +u(n s + 1/2) (18) 

where w — Co(\p v \/pFr)^/2fi/\fi x \. One can see that the 
Eq. dTSl) yields e > \fi x \. Hence the energy branches at 
fi x < can have arbitrary high energy when fi x — > —00. 
This tendency can be seen in the Fig. (J^b) ■ 

Besides the energy branches determined by Bohr- 
Sommerfeld quantization with n s > 1 there exists a so- 
called 'zeroth branch' of the spectrum which contains the 
zero modes e — 0. This branch can not be obtained from 
the semiclassical approach therefore the different treat- 
ment is needed^. 



2. Quasiclassical approximation. 

To find the zeroth branch analytically we assume the 
restriction (p x /pf)Px£, 3> 1 where p x — y/2m,fi x is a con- 
stant number. Qualitatively this restriction corresponds 
to the regime (iii) in Fig|3]when the deviation p x (x) of 
momentum p x (x) = \J1mfi x + p x {x) along the orbit is 
much smaller than its average value \p x | <C \J1mp, x . In 
this case we can neglect the second order terms in p x still 
treating it as an operator so that [e{p) — fi] ps —iv x d x 
where the projection of the Fermi velocity on the x axis 

is v x = p x /m = ±^Jp F — pi — p y /m. This allows to 

consider the regimes when the closed orbits become so 
small that semiclassical approximation does not apply. In 
particular it allows to calculate the zeroth branch which 
appears when the distance between Andreev reflection 
points is much smaller than the coherence length £. In 
this way one obtains a system of abridged equations for 
the wave function envelopes which is called the system 
of quasiclassical Andreev equations. The quasicalssical 
approximation inherently misses the normal reflection of 
particles from the order parameter inhomogeneities. As 
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FIG. 4: (Color online) (a) Closed orbits p x = p x (x) deter- 
mined by the semiclassical Eq. (|14p at the domain wall for 
p z = 0, p y = 0.5pf and co/vf = 2.5 10~ 3 . Shown by 
red dashed, green solid and blue dash-dotted curves are the 
regimes (i), (ii) and (iii) discussed in the text. The posi- 
tions of Andreev x a i,2 and normal x n i,2 reflection points on 
the semiclassical orbits is indicated, (b) Comparison of energy 
spectrum given by quasiclassical and semiclassical approxima- 
tions. The quasiclassical spectrum given by Eq. ()24p is shown 
for n q = 0, 1 by dashed red and blue lines. Quasiclassical ap- 
proximation takes into account only Andreev reflection hence 
the spectrum is doubly degenerate by sign(j> x ). The semiclas- 
sical approach describes both Andreev and normal reflection 
hence the degeneracy by sign(p x ) is removed. The branches 
given by Bohr-Sommerfeld quantization (|17|) of semiclassical 
orbits with n s — 1,2 are shown by solid magenta lines. The 
number of semiclassical branches is doubled due to the nor- 
mal reflection when the Andreev orbits shown in the panel 
(a) by green solid lines merge into the single one shown by 
blue dash-dotted line with both normal and Andreev reflec- 
tion points. 

we will see below the normal reflection is crucial to de- 
scribe the Fermi arcs and Majorana states at the domain 
wall considered. However away from the Fermi arcs the 
quasicalssical approximation matches the exact spectrum 
obtained numerically. 

To obtain analytical solutions of Andreev equations 
first let us introduce the transformation of the quasipar- 
ticle wave function tp — (it, v) T 

g+ = v + u (19) 
g- = v - u (20) 

so that the system of Andreev equations acquires the 
form 

( H+ £ \ CoP A( 9+ )=0 (21) 
\e-cop x h- J \ g _J 

where h± = —iv x d x ± icof(x)p y . 

The system (f2~Tj) can be transformed to the decoupled 
second-order equations 

[dl + U 0+ cosh" 2 x] g+ = Eg+ (22) 
[d 2 x + U - cosh" 2 x] g- = Eg- (23) 

where E = (r/£) 2 + a 2 - (p F e/p x A ) 2 , U 0+ = a + a 2 , 
a = rc p y /v x and U -(p y ) = U 0+ (—p y ). The Eq.(g2]) 
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has eigenvalues 

E nq = (1 + 2n q - Vl + 4C/o+) 2 /4 - (n, - a) 2 

where n q > is integer, which results in the quasiclassical 
spectral branche o 17 i 25 

e„ 9 - ± A M ./2cm, - n 2 + (r/0 2 - (24) 

Note that the energy branches (|24p are twofold degen- 
erate with respect to the sign change of p x projection. 
The zeroth spectral branch is given by the Eq.(|24|) with 

n q = 

£o(Py,Pz) = c sign(p y )\p x \. (25) 

The dependencies of zeroth branch (f2"5"j) £o = £o{Hx) at 
fixed p y and eo = £o(p y ) at fixed p z are shown in Figj5] 
by red dashed lines and in Fig|6] by blue dashed lines 
correspondingly. At p y — this spectral branch is dis- 
continuous and merges the edge of continuum e c = ±cqp* 
at p y — ±0. 

Substituting the energy of zeroth branch into the 
Eq. ([T6l) for the coordinates of Andreev reflection points 
we obtain \x a i :2 \ = Pp 1 ^ /£)(pF /\p v \) ~ Pp 1 provided 
\p y \ is not too small. For such a small distance between 
reflection points the semiclassical approximation is not 
valid and therefore it can not describe the zeroth branch. 

The spectral branches with n q > given by Eq. (f24"]) 
exist only in the limited range of parameters a/p_f/£ <C 
\p x \ < P* x where p% = \p y \n q /[n 2 q - (r/£) 2 ] provided that 
n q > r/£. At p x = p* the branch merges continuum. 
Thus behaviour us demonstrated in FigsHfc where the 
energy branches (|24|) with n q = 0, 1 are shown by dashed 
red and blue lines as function of (j, x = /i — (p 2 + p 2 )/2m 
for the parameters r = £/2 (width of the domain wall), 
cq/vf = 2.5 10~ 2 and p y = 0.5pf- The zeroth branch 
merges continuum at p y = where the spectrum Eq. (|25p 
is discontinuous. This behaviour is shown by dashed blue 
lines in Fig.©. 

3. Exact spectrum near the Fermi arcs. 

As pointed out by Nakahara 25 all quasiclassical 
branches (1231) formally contain zero energy states at 
p x = but Eq.(f2"5"j) is applicable if (Px/pf)Px£ ^S> 1. In 
the limit p x — > and arbitrary values of p y and p z the 
spectrum can be found numerically by solving the eigen- 
value problem for the Eq. ([21]) . The resulting several 
lowest energy spectrum branches e = e(/x x ) are shown in 
the Figj5] by solid red and blue lines as the function of 
Hx = A* — (Py + p 2 z )/2m. Remarkably the twofold degen- 
eracy p x — > —Px of the quasiclassical spectrum is broken 
at < n x <C fi and at fi x < 0. Therefore each of the qua- 
siclassical energy branches (|25|) found by Nakahara splits 
by two modes as fi x +0 (and p x — \/2mfi x — > 0). 
The splitting is caused by normal reflection of quasipar- 
ticles from the order parameter spatial inhomogeneity. 



The normal reflection is particularly important for the 
quasiparticle trajectories passing almost parallel to the 
domain wall plane when p x — > 0. In this case the nor- 
mal reflection leads to the drastic change in the spectral 
branches behaviour near the zero energy point shown in 
Fig E] (compare the dashed and solid lines). The influence 
of normal reflection on the spectrum of bond fermions 
was analyzed above in the framework of semiclassical ap- 
proach in the section pil B Comparing the Figs.(]4jD) 
and ((5|) one can see that the semiclassical approach qual- 
itatively describes the doubling of spectral branches due 
to the normal reflection but gives rather large quantita- 
tive discrepancy with the exact spectrum. 

In accordance with general topological argument in 
Sec. ()III A[) the only one zeroth branch remains at each 
fixed value of momentum projection p y (and spin projec- 
tion a x for the domain wall of the A phase) . Exact zeroth 
branches for p y = ±0.5pF are shown in the Figj5]by thick 
red solid lines. The spectrum in Fig(5] is invariant with 
respect to the transformation e(p y ,a) = —e(—p yi —o) 
where a = a x for the case of domain wall and a = <r± 
for that of the AB interfere. Note that in contrast to 
the quasiclassical zeroth branch shown by red dashed 
line the exact zeroth branch intersects at the Fermi level 
at finite value of \i x = fi x0 > 0. This crossing point 
determines the positions of Fermi arcs shown in the 
Fig. (Eb) by blue and green solid lines. For example the 
Fermi arcs cross the line p z = at p y — ip y a where 

PyO = y/2m(lJ, - Uxo) < PF- 

Near the projection of Weyl point on p y p z plane in 
the limit p y — > and \p z \ — » pf the behaviour of zeroth 
branch can be found analytically using the approach of 
Ref . (|20T ) . We will treat the first term in the Hamiltonian 
(fljl as perturbation. The rest of the terms form zero order 
Hamiltonian 

H = tiPx - T 2 p y f(x) (26) 

where p x — —id x . We assume the model form of the do- 
main wall f(x) = tanh(x/r). The Hamiltonian ([26| has 
zero energy eigen state with the wave function compo- 
nents ?/>o = {uq,v ) 

u = (27) 
v = N- 1/2 cosh- a (x/r) (28) 

where a — rp y and N = J X Qo cosh^ 2a (x/r)dx. The per- 
turbation of the zero energy level is given by 

f 

£o(Px,Py) = ^-Vx (29) 
where p 2 = (iPq\p x \ipo)- The Fermi arc is then given by 

2 2 -2 2 /on\ 

Pz = Pf ~ Px ~ P y (30) 

In the limit p y — > we obtain p 2 — p 2 so that the Fermi 
arc is given by p 2 — p 2 F — 2p 2 and ends at p z = ±pp 
in accordance with general topological arguments. The 
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Eq. (l29l) describes analytically the zeroth branch of the 
spectrum shown by thick red lines in Fig. ([3]) near the 
projections of Weyl points when p y — > 0. 

The splitting of quasiclassical zeroth branches e = 
£o(Py) at fixed value of p z — is shown in Fig [SI Red 
solid lines is the exact spectrum obtained numerically and 
blue dashed lines are the quasiclassical zeroth branches 
(f!?5)) . Away from the level e = the correspondence of 
quasiclassical and exact spectra is of a very good accu- 
racy. The splitting takes place in the small vicinity of 
p y = ±p* where the quasiclassical branch intersect the 
Fermi level. The behaviour of spectral branches near 
this point is shown in the zoom inset in the Fig.®. All 
branches are discontinuous at p y = where they merge 
the continuum. Shown in Fig. ([6]) the spectrum of bound 
fermionic states at the domain wall in 3 He-A at p z = 
coincides with the spectrum of electronic states at do- 
main walls in p + ip superconductor SV2-RUO4. 




FIG. 5: (Color online) The spectrum of surface states at the 
A phase domain wall and the AB interface e = e(jj, x ) where 
fJ-x = fJ- — (Py + pi) /2m. The parameters are r = £/2 (width 
of the domain wall) and cq/vf = 2.5 1CP 2 , \p v \ = 0.5pf- 
Shown by solid red and blue lines are the spectral branches 
determined by the exact solution of BdG Eq.{TJ. Thick red 
solid lines show the zeroth branches at p y — ±0.5pF- Dashed 
lines represent quasiclassical energy branches (I24|l for n q — 
(red line) and n q — 1 (blue line) degenerate by sign(p x ). The 
degeneracy by sign(p x ) is removed at the region \p x \ <C ft 
and fi x < where each quasiclassical branch splits by two 
modes. In accordance with general topological argument the 
only one zeroth branch remains at fixed value of p y (and spin 
projection a x for the domain wall of the A phase) shown by 
the red solid line. The spectrum of delocalized states is shown 
by yellow shading and the edge of continuum by thin dotted 
line. 
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FIG. 6: (Color online) The lowest energy spectral branches 
of the surface states at the domain wall of the A phase and 
at the AB interface as functions of p y at fixed p z . Shown 
by blue dashed line is the quasiclassical zero energy branch 
(|25[) which ends at p y = ±p y where p y = \/p 2 F — pi- Exact 
branches shown by solid red lines continue at \p y \ > p y . At 
p y = the discrete spectrum merges with continuum shown 
by the yellow shading. Inset: zoomed area demonstrating in 
detail the splitting of zeroth branches and the behavior of 
exact spectral branches near the end point of quasiclassical 
spectrum branch. 



IV. DISCUSSION 

In contrast to the Fermi arcs on 3 He-A surface^ the 
Fermi arcs on the domain wall in 3 He-A and on the AB 
interface can not be obtained in the framework of qua- 
siclassical approximation. Indeed the quasiclassical ap- 
proximation when applied to these systems yields that 
any of the subgap spectral branches given by Eq. ([24|) 
intersects the Fermi levels. Thus the number of zero 
modes in quasiclassics is model depended and does not 
satisfy the index theorem. The reason behind this dis- 
crepancy is that the zero modes exist for p x = (see 
dashed lines in FigfJJ}) which corresponds to the quasi- 
classical trajectories almost parallel to the plane of do- 
main wall when one needs to get into account the normal 
reflection of qasiparticles from the order parameter inho- 
mogeneity along the x axis. Since the normal reflection 
is inherently missing in quasiclassical Andreev equations 
(without diagonal potential) this approximation fails to 
describe zero modes and Fermi arcs on the domain wall 
in 3 He-A and the AB interface. 

The modification of quasiclassical results due to the 
normal reflection can be analyzed by employing the semi- 
classical approximation (see Section UlI B II) which yields 
three possible types of closed orbits in (p x , x) phase space. 
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In case when quasiclassical approximation is valid the 
regime (ii) of only Andrcev reflection is realized with 
two separate orbits at p x > and p x < shown by the 
green solid lines in Fig. 0^. The Bohr-Sommerfeld quan- 
tization (|17p yields the same energy of bound states for 
the two orbits which is degenerate by sign(p x ) and coin- 
cides well with quasiclassical result for n q > 0. However 
for smaller fj, x two Andreev orbits merge into a single 
one shown in Fig|4ji by blue dash-dotted line. This orbit 
has both Andrcev and normal reflection points and there- 
fore is not described by the quasiclassical approximation. 
Moreover merging of Andreev orbits doubles the orbit 
area. Therefore the number of subgap states accord- 
ing to Bohr-Sommerfeld quantization is doubled which is 
manifested by the appearance of the second semiclassi- 
cal branch in FigfJb simultaneously with merging of two 
Andreev orbits in Fig|4^,. In exact spectrum obtained 
by solving numerically the BdG system the branch dou- 
bling is manifested in splitting of quasiclassical branches 
in Figs l5l6l Thus the branch splitting is the direct re- 
sult of normal reflection of quasiparticles moving almost 
parallel to the domain wall. 

Affected by the enhanced normal reflection the be- 
haviour of the energy branches changes in p x — > limit 
qualitatively compared to the quasiclassical result. At 
fixed value of momentum projection p y (and spin projec- 
tion a x in case of A phase domain wall) there is only one 
branch left which intersects the Fermi level and supports 
Majorana states in accordance with the topological in- 
dex theorem. The rest of spectral branches become non- 
monotonic at small values of \p x \ = ^2mfi x at /i x > 
and turn upwards at fi x < as can be seen in the Fig{5j 

V. CONCLUSION 

We demonstrated that Fermi arcs on the domain wall 
in the A-phase and at the interface between the A-phase 
and the B-phase of superfluid 3 He obey the index theo- 
rem, which connects topological properties of bulk states 
on two sides of the walls with topology of zero energy 
bound states. In contrast to the other systems support- 
ing Fermi arcs the bulk-edge correspondence here can be 



established only beyond the quasiclassical approximation 
by solving the exact BdG system numerically. The do- 
main wall in A-phase contains 4 Fermi arcs, if one takes 
into account spin degrees of freedom, while the AB in- 
terface contains 2 Fermi arcs, i.e. the same number as 
the surface of the A-phase. The reason for this difference 
is that the Fermi arc is determined by the Weyl points 
in the bulk states on two sides of the wall. The 3 He-B, 
though being a topological superfluid, does not contain 
Weyl points and thus does not contribute to the number 
of the Fermi arcs. Note that the spectrum of localized 
fcrmions at the domain wall of the A phase is identical 
to the spectrum of electronic states bound at the domain 
wall in p + ip superconductor Sr 2 Ru04. 

Majorana fermions living on the Fermi arc at the AB 
interface may give an additional contribution to the cal- 
culated friction force acting on the moving interface 3 - 3 - at 
very low temperatures. This is important for the devel- 
opment of the Kelvin-Hclmholtz instability of the moving 
AB interface observed in Ref. [34] (see also Chapter 27 of 
Ref. d). The latter instability has an analogue with the 
ergoregion instability of black holes discussed in Chapter 
32 of Ref. 0. 

We considered the simplest most symmetric realiza- 
tions of the A- A and A-B interfaces. The problem for 
future investigations is weather the Fermi arc survives 
or not, if the symmetry of the interface is violated. For 
that the relative homotopy group formalism applied to 
Green's function 3 ^ is probably required. 
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